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of global Z3 symmetry, and possible splitting of deconfinement and chiral transitions in a strong 
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I. INTRODUCTION 

The phase diagram of QCD has a very rich structure, 
and keeps unveiling new, sometimes unexpected phases 
of strong interactions as one plays with thermodynamic 
control parameters such as the temperature, chemical 
potentials and quark masses (for recent results, see e.g. 
Ref. [1]). A particular external control parameter that 
is often overlooked, however, is the magnetic field. Since 
quarks and charged hadrons couple to the magnetic field 
and neutral hadrons do not, this control parameter can 
affect the phase structure of QCD in a nontrivial fash- 
ion. Moreover, it is present in most physical systems 
exhibiting deconfinement and chiral symmetry restora- 
tion. The example that was considered to be the most 
spectacular is from astrophysics, provided by magnetars 
[2]. Besides, magnetic fields are expected to be strong 
and play an important role in structure formation in the 
early universe, during the epochs of the electroweak and 
the QCD primordial phase transitions [3]. However, it 
was recently found that non-central high-energy heavy 
ion collisions might generate the most intense observed 
magnetic fields, much stronger than in magnetars, reach- 
ing values of _B ~ 10^^ Gauss, which corresponds to 
eB ~ 6to^, where e is the fundamental charge and is 
the pion mass. These magnetic fields are short-lived for 
very high energies but play an important role in possible 
experimental signatures of strong CP violation and the 
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phenomenon of the chiral magnetic effect [JHZ]- 

In this paper we present a general picture for the 
temperature-magnetic field phase diagram of QCD, re- 
stricting our analysis to a vanishing baryonic chemical 
potential. In particular, we investigate how the chiral and 
the deconfining transitions are affected, and uncover new 
phenomena such as the paramagnetically-induced break- 
ing of Z3. We also demonstrate the possible splitting of 
the deconfinement and chiral transition in an external 
magnetic field. 

For this purpose, we adopt the linear sigma model cou- 
pled to quarks with two fiavors, Nf = 2, and to Polyakov 
loops (PLSMq) under the effect of a magnetic background 
field as an effective theory for the thermodynamics of 
QCD. The linear sigma model coupled to quarks ^ has 
been widely used to describe different aspects of the chi- 
ral transition, such as thermodynamic properties and the 
nonequilibrium phase conversion process [97:20j . It has 
also been combined to Polyakov loops to include confine- 
ment [2TI - [24] . However, as will be clear in the next sec- 
tion, we perform this generalization in a different fashion, 
which is also more natural in the presence of an external 
magnetic field. 

In our effective model, PLSMg, we compute and dis- 
cuss each one-loop contribution to the effective poten- 
tial for the approximate order parameters of the chiral 
and deconfining transitions, given by the chiral conden- 
sate and the expectation value of the Polyakov loop, re- 
spectively. We treat the external magnetic field, that 
plays the role of a thermodynamic control parameter, 
as a constant and uniform field, and make no assump- 
tion on its intensity. In Fig. [l] we show a cartoon of the 
temperature-magnetic field phase diagram one would ex- 
pect from previous results for the limits of strong and 
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weak fields, for chiral and deconfining transitions treated 
separately, as discussed below. 

Expected phase diagram 



Restored chiral symmetry, / 
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Magnetic field, eB 

FIG. 1. Expected magnetic field-temperature phase dia- 
gram of QCD. Tiie thick lines indicate first-order transitions, 
the filled circles are the (second-order) endpoints of these 
lines, and the thin dashed lines stand for the corresponding 
crossovers. A new phase with broken chiral symmetry and 
deconfinement appears at high magnetic fields. 

Modifications in the vacuum of strong interactions 
by the presence of an external magnetic field have 
been investigated previously within different frameworks, 
mainly using effective models [25M33] , especially the NJL 
model [33], and chiral perturbation theory [55H37) . but 
also resorting to the quark model |38j and certain lim- 
its of QCD |39j. Most treatments have been concerned 
with vacuum modifications by the magnetic field, though 
medium effects were considered in a few cases, as e.g. in 
the study of the stability of quark droplets under the 
influence of a magnetic field at finite density and zero 
temperature, with nontrivial effects on the order of the 
chiral transition (JHI- Magnetic effects on the dynamical 
quark mass [U] , as well as magnetized chiral condensates 
in a holographic description of chiral symmetry breaking 
[32] J were also considered. Lattice simulations were also 
used to demonstrate an enhancement of the chiral con- 
densate and the appearance of chiral magnetization of the 
chirally broken vacuum under the influence of a strong 
magnetic field J43J . Ring-diagram corrections to the ef- 
fective potential in a weak magnetic field were considered 
in Ref. [331 . Recent applications of quark matter under 
strong magnetic fields to the physics of magnetars using 
the NJL model can be found in Ref. |49i 150j . 

The case of the thermal quark-hadron transition in a 
magnetic field was studied recently in Ref. [51]. It was 
found that the critical temperature and the latent heat 
are both diminished in the presence of the magnetic field, 
and that the first-order deconfining turns into a crossover 



at sufficiently strong magnetic fields. On the other hand, 
it has been shown in Ref. [521 ^^^^ chiral transi- 
tion can be dramatically modified in the presence of a 
strong magnetic field. In particular, the smooth crossover 
found for vanishing baryonic densities can be turned into 
a first-order phase transition if the magnetic field is suf- 
ficiently high, the threshold being close to fields that can 
be achieved in current and future heavy ion collision ex- 
periments. For the chiral transition - on the contrary to 
the case of the deconfining transition of Ref. [51] - the 
critical temperature seems to increase with the magnetic 
field, thus favoring a first-order behavior, which goes in 
line with the theoretical expectation on the role played 
by magnetic fields in the nature of phase transitions [53] . 
On the other hand, at very weak magnetic field both 
confinement and the chiral transitions should again be 
crossovers [54] . This justifies our preliminary approxi- 
mate cartoon of Fig. [l] 



More recently, the chiral magnetic effect was also stud- 
ied in the context of the PNJL model [6 . Here we pro- 
pose a different effective theory for the study of the phase 
structure of QCD in the presence of an external magnetic 
field, the linear sigma model coupled to quarks and to the 
Polyakov loop, and compute the T—eB phase diagram in 
the presence and, separately, in the absence of vacuum 
corrections. The contribution from these corrections is 
very large, though the issue of their relevance or neces- 
sity is non-trivial in this context as also in different ones, 
such as the Higgs potential in the electroweak theory and 
the Walecka model in high-density nuclear physics |55| . 
We choose to present results in both scenarios and ad- 
dress the issue in the end. Nevertheless, we believe that 
lattice results will be crucial in setting the most appro- 
priate effective theory approach. The PLSMg can, of 
course, also be used to study chiral symmetry breaking 
and deconfinement in the absence of a magnetic field. 
We present a few results in this case but, keeping our 
focus on the phase structure in the presence of a mag- 
netic background, leave a thorough study for a future 
publication ^56^. 



The paper is organized as follows. In Section [IT] we 
provide a description of the effective theory that allows 
us to study the influence of an external magnetic field 
on confining and chiral properties of QCD, the linear 
sigma model coupled to quarks and to the Polyakov loop 
(PLSMq). We discuss the quark, chiral and confining sec- 
tors, the total free energy and the connection to various 
observables. In Section [III] we compute all one-loop con- 
tributions to the effective potential, discussing vacuum 
and thermal corrections in detail. Section IIVI contains 
our results for the effective potential and a discussion of 
the temperature-magnetic field phase diagram. Section 
Ivl contains our conclusions and outlook. 
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II. EFFECTIVE THEORY 



A. Quark, chiral and confining variables 



of the external magnetic field on color confinement and 
chiral symmetry breaking simultaneously. 

We can write the Lagrangian of PLSM^ as follows: 



The effective field theory that we denote by PLSM^ 
has three types of variables: the quark field ip{x), the 
chiral field 4>{x) and the Polyakov loop variable L{x). 

The confining properties of QCD are described by the 
complex- valued Polyakov loop variable L, which plays 
the role of an approximate order parameter for the con- 
fining phase transition in QCD. The Polyakov loop is a 
singlet part of the 3x3 matrix $ ("untraced Polyakov 
loop") which belongs to the adjoint representation of the 
SU{3) color group, i.e.: 



$ = 7? exp 



l/T 

i J dTA4ix,T) ,(1) 




where — iAo is the matrix-valued temporal compo- 
nent of the Euclidean gauge field and the symbol V 
denotes path ordering. The integration takes place over 
compactified imaginary time r. 

The expectation value of the Polyakov loop L is an 
exact order parameter of the color confinement in the 
limit of infinitely massive quarks: 



Confinement 



(L) = 



low T 
high T 



(2) 



The chiral features of the model are encoded in the 
dynamics of the 0(4) chiral field 



= (o-, tt) , 



(3) 



where tt is the isotriplet of the pseudoscalar pion fields 
and a is the chiral scalar field which plays the role of 
an approximate order parameter of the chiral transition. 
The expectation value of the field a is an exact order 
parameter in the chiral limit, in which quarks and pions 
are massless degrees of freedom: 



Chiral symmetry 



(a) 



+ 




low T 
high T 



(4) 



The up and down fields of the constituent quarks are 
grouped together into the doublet fcrmion field 



(5) 



which plays a central role in our discussion. The fermion 
field couples two other variables, the Polyakov loop L 
and the chiral field to each other, thus linking confin- 
ing and chiral properties together. The quarks ^ are also 
coupled to the external magnetic field since the u and d 
quarks are electrically charged. Since the quarks are also 
coupled to the Polyakov loop and the chiral field ^, it is 
clear that the external magnetic field will affect the chi- 
ral dynamics and the confining properties of the model. 
Thus, this model allows us to investigate the influence 



L = Lq{i,,i,,4>,L) + L^{<i>)^ Ll{L). (6) 

B. Quark sector 

The first part of the full Lagrangian ([6]), 

= -0 [i^ - g{o + il-^T • 7?)] -0 , (7) 

describes the constituent quarks, which interact with the 
meson fields 0, the Abelian gauge field and the SU{3) 
gauge field via the covariant derivative: 



(8) 



The Abelian gauge field describes the influence of the 
external magnetic fleld B aligned along the third direc- 
tion, Bi — B5iz, for convenience: 



(a°,a) = (0,-Sy,0,0), 



(9) 



and the electric charges of the quarks are deflned by the 
following matrix: 



_ ( qu 
qa 



-I' 




(10) 



where e is the elementary electron charge. 

The SU{3) gauge field A^ represents a nontrivial back- 
ground due to the Polyakov loop ([T]). It is convenient to 
diagonalize the untraced Polyakov loop $. In this gauge 
the field Ai, which enters the quark Lagrangian ([t]), is 
diagonal: 



Aa 



h ^4 



(3) 



UA 



(8) 
4 J 



(11) 



where fs.s = A3_8/2 are the Cartan generators of the 
SU{?>) gauge group. The S'C/(3) gauge field is prescribed 
to take constant values, being linked to the Polyakov loop 
as follows: 



$ = exp 



diag (e'^^S e*^% e^'^^) , 



(8)- 



(12) 



with 



V3 



AT 

T 
Af 

T 
Af 

T 



_ in 4(3) I 1 a(^)\ 

- T\2^i ^ 2^3 4 ) : 

- T\ 2^4 + 273^4 J 
_ _1J_4(8) 

^ T v/3 4 ■ 



(13) 
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C. Chiral Lagrangian 

The chiral part of the PLSM, Lagrangian is given by 
the second term of Eq. that can be written in the 
Minkowski space-time as foUows: 



(14) 



Here we have introduced the charged tt* and neutral 
mesons, respectively, 



v2 



(15) 



The electric charges of the 7r*-mesons are ±e, so that the 
covariant derivative in Eq. ( 14 ) reads 



(16) 



This derivative does not involve the Polyakov loop, con- 
trary to the covariant derivative acting on quarks ([s]), 
since the pions are colorless states. 

In Ec^. ( 14 ) , stands for the potential of the chiral 
fields. This potential exhibits both spontaneous and ex- 
plicit breaking of chiral symmetry: 



(17) 



2 2 



In the second line of this equation we explicitly show 
the quadratic form of the potential which provides the 
masses to the meson fields. Here, the masses mo- and 
correspond to the vacuum masses of the sigma and 
the pion mesons, which are used to fix the parameters A 
and V in the classical potential, since /^r is also known, 
as explained below. 

As customary in the linear sigma model, we follow a 
mean-field analysis in which the mesonic sector is treated 
classically whereas quarks represent fast degrees of free- 
dom. The parameters of the Lagrangian ( 14 ) are chosen 



to match the low-energy phenomenology of mesons in the 
absence of magnetic fields and at zero temperature, i.e. 
in the vacuum [57] . Using the notation of Ref . [T^ , this 
implies the following expectation values for the conden- 
sates: (cr) = and (tt) = 0, with « 93MeV and 
« 138 MeV. The masses in Eq. (17) should be con- 



sidered as parameters of the model, which coincide with 
corresponding meson masses at T = and B = 0. 

The explicit symmetry breaking term, as usual, is de- 
termined by the PC AC relation h — fTrin^, so that 
v'^ = — m^/X and = 2A/^ -I- in the vacuum. 
Choosing the value of the quartic interaction A = 20, 
one gets a reasonable cr- mass ~ 600 MeV. The con- 
stituent quark mass is given by ruq = mq{{a)) — g{<y), 
and, choosing g — 3.3 at T = 0, one obtains for the 
constituent quarks in the vacuum sa 310 MeV. At 
low temperatures the quarks are not excited, and the 



model ( |14| reproduces results from the usual linear cr- 
model without quarks [55j . 

If the explicit symmetry breaking term is absent, hq — 



0, then the model (14) experiences a second-order phase 



transition 15.8] from the broken symmetry phase to the 
restored phase at a critical temperature Tc = y/2v. The 
presence of the explicit symmetry breaking term changes 
the phase transition into a smooth crossover. 

The pion directions play no major role in the process 
of phase conversion we have in mind, as was argued in 
Ref. [13] for the case of the chiral transition. Therefore 
in what follows we focus on the sigma direction of the 
chiral sector. However, the coupling of charged pions 
to the magnetic fields might be quantitatively important 
in a possible development of charged chiral condensates. 
This issue, though, makes the computation technically 
much more involved and will be addressed in a future 
publication |56j . 



D. Confining potential 

We choose to couple the Polyakov loop to the linear 
sigma model in the same fashion as was successfully im- 
plemented in the case of the Nambu-Jona-Lasinio model 
in Rcfs. [55ti61| . The confining properties of the system in 
our model are accounted for by the third term in Eq. (|6| , 
which describes the Polyakov loop variable ([!]): 



-Vl{L,T). 



(18) 



We treat the Polyakov loop in a mean-field approach, 
neglecting a kinetic term for this variable. Although the 
Polyakov loop is, strictly speaking, a static quantity, a 
kinetic term is often added to describe the dynamics in 
an effective model for deconfinement [BH [53] . 

Notice that, in our approach, the Polyakov loop vari- 
able L is not coupled to the chiral variables (f) in the bare 



Lagrangian (18). The reason is that the chiral variables 



are colorless degrees of freedom which should not affect 
the Polyakov loop at least in a first approximation. Of 
course, the mesons (f> are made of colored quarks which 
do feel colored degrees of freedom like L. Thus, the in- 
teraction between the chiral and confining variables will 
inevitably appear later, after the quarks ip are integrated 
out. This approach is not unique. In fact, it is com- 
mon to couple the Polyakov loop directly to chiral fields 
phenomenologically as a way to link the chiral and de- 
confining transitions in a generalized Landau-Ginzburg 
theory, in which the behavior of an order parameter in- 
duces a change in the behavior of non-order parameters 
at the transition via the presence of a possible coupling 
between the fields [2Tti23l. 



The term ( 18 ) is essentially the potential term of the 



Polyakov action in the pure Yang-Mills theory not cou- 
pled to quarks. Quark effects manifest via coupling with 
the quark matter field in the first part of our Lagrangian, 
Eq. (0. 
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ties 



The effective potential (18) has the following proper- 
it should satisfy the center Z3 symmetry, 



btiL 



-1.0 



-0.5 



0.0 



0.5 



1.0 



0,1,2. 



(19) 



which is realized in the limit of the pure gauge theory 
(the coupling to dynamical quarks break this symmetry). 
According to Eq. the potential (18) should have an 
absolute minimum at L = in the confined phase. Re- 
sults from lattice simulations [M] show that the decon- 
fining transition in SU{5) gauge theory takes place at a 



T = 0.8 To 



-1.0 



-0.5 



critical temperature T< 



SU{3) 



270 MeV. As the system 



0.0 

ReL 



0.5 



1.0 



undergoes the phase transition the single minimum at 
L = of the potential splits into three degenerate min- 
ima labeled by the Z3 variable. As a consequence, in 
the deconfined phase the Z3 symmetry is spontaneously 
broken, and (L) ^ 0. 

Following Ref. [ST] we use a specific form for the phe- 
nomenological potential of the Polyakov loop: 



IiiiL 



VLiL,T) 



--a{T)L*L 



rp4 

-b{T) In 1-6 L*L + 4 ( L*^ + ) - 3 {L*Ly 



(20) 




with 



a{T) = ao + fli 



02 



6(T) 



(21) 
(22) 



FIG. 2. The phenomenological potential of the Polyakov 
loop ((20} at T = 0.8 To (top) and T = 1.2 To (bottom). 



E. Total free energy and observables 



where Tn is the critical temperature in the pure gauge r r , 1 • ■ i, 

^ f b b j^j^g j-j-gg gngj-rry ot the system is given by 

case: Tq = Tsuis) = 270 MeV. B.y j & j 



The phenomenological parameters in Eq. ( 22 ) are 



ao 

0.2 



16 7rV45 : 
15.2, 



3.51, 



ai = 
63 = 



—2 47 

-175 : (23) 



T 

n = -—inz, 



(24) 



Consistently with the original definition ([ij, the poten- 
tial (20) limits the values of the Polyakov loop to the 



where V^d is the spatial volume and Z is the partition 
function 

i/r 



interval L* L ^ 1 because of a logarithmic divergence in 
(20). The value L*L = 1 is reached in the limit T — cx). 



VipVipVaVn exp 



In the confined phase, T < Tq, the potential has one triv- 
ial minimum [as illustrated in Figure [2] (top)], whereas in 
the deconfined phase, T > Tq, the potential has three 
degenerate minima [see Fig. [2] (bottom)]. The latter are 
mutually connected by Z3 transformations. 

The set of parameters ([23]) is obtained by demanding 
the following requirements valid for pure SU (3) gauge 
theory [61 : 

(i) the Stefan-Boltzmann limit is reached at T — > 00; 

(ii) a first-order phase transition takes place at T = Tq; 

(iii) the potential describes well lattice data for the 
Polyakov loop and for the thermodynamic func- 
tions such as pressure, energy density and entropy. 



-i J dt J (TxC 




(25) 



In the mean field approximation, treating the meson 
fields classically and keeping only the sigma direction in 
the mesonic sector, the free energy is given by 

n[a, L- T, B) = V^{a, tt ) + Vl (L, T) + n^ia, L, T) , 

(26) 

where the potential for the chiral field is given in 
Eq. (17) and the potential for the Polyakov loop is 



defined by Eq. ( 20 ). In ( 26 ) we explicitly show the depen- 



The uncertainties for the parameters in Eq. (23) are: 



6 % for fli , less than 3 % for a2 , and about 2 % for 63 [5T] . 



dence of the free energy on the external magnetic field 
the temperature T, and the mean field value of the chiral 
field (7, as well as the expectation value of the Polyakov 
loop L. 

The temperature and the magnetic field impose the 
fixed external conditions whereas the expectation values 
of the dynamical fields are obtained by minimization of 
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the free energy ( 26 1 with respect to the dynamical vari- 
ables: 



dX 



= 0, 



X = tj, L. 



(27) 



The mass of the cr-mesons is given by the curvature 
of ri evaluated at the global minimum (27) of the free 
energy, i.e.: 



(28) 



Notice that the mass (281 depends implicitly on the 



temperature T and the strength of the magnetic field 
B. The zero-temperature and zero-field value of the a- 
meson mass (defined in the preceding Section) is related 
to Eq. (pSl) as follows m„ = m„{T = 0, B = 0). 



III. FREE ENERGY AT ONE LOOP 

We compute the free energy of the system in the one- 
loop approximation. To this end, we functionally inte- 
grate over the quark fields and take the integral over 
the fluctuations of the mesonic degrees of freedom near 
the minimum of the free energy, determined by Eq. (27 1. 



These integrals yield thermal determinants, either in the 
background of the magnetic field B or for vanishing field. 

The one-loop correction coming from quarks can be 
written in the form: 



exp 



V^pVip exp 



n 



detT [ip^''^ - g{a + ij^r ■ n)] 
det [ip - g{a + ij^r ■ tt)] 



(29) 



where the quark Lagrangian Cq is given in Eq. (0. The 
notation "det^" means that the determinant is taken at 
nonzero temperature, while "det" indicates that we con- 
sider the T = case. 

The numerator in Eq. ( 29 ) involves the covariant 
derivative ([s]) because we calculate the contribution of 
the fermion loop in the presence of an external magnetic 
field and in the background of the nontrivial Polyakov 
loop. The magnetic field B and the Polyakov loop L 
affect the chiral degrees of freedom via the covariant 
derivative. Thus, the integration over quark fields pro- 
vides an effective potential for the meson fields, a and 
TT, and the Polyakov loop, L. The coefficients of this po- 
tential should in general depend on the strength of the 
magnetic field B and on the temperature T. The deter- 
minant in the numerator of Eq. ( 29 ) is regularized by a 



similar determinant in the denominator which is calcu- 
lated at T = 0, S = and in the absence of the confining 
background. Thus, the quark one-loop correction is ex- 
pressed in terms of a ratio of fermionic determinants. In 
what follows we assume that the vacuum of the model is 
defined by (cr) = /tt, (tt) = 0. 



Following Ref. [5T] it is convenient to regroup in 
Eq. (29) the ratio of the determinants, and represent the 



free energy of each species of the light quarks as follows: 



det(i^(9) - m,) 



det(zj9 — TUq) 

detT(j^^'^ - mq) 
det(i^('?) - fUq) 



(30) 



where q stands for either u or d. The first ratio corre- 
sponds to the response of the quark loops to the external 
fields at zero temperature. This term represents the vac- 
uum contribution, 51™'^. The second ratio in Eq. (30) 



represents the finite-temperature correction due to ther- 
mal fermion excitations. As in Ref. [5T], we call this 
contribution the paramagnetic part of free energy, f^P''''^. 

In this paper, as mentioned previously, we consider 
two approaches for the PLSMq, one that includes in- 
dependent vacuum effects and one that does not. We 
intentionally ignore pure vacuum corrections, in order to 
compare our results with what is found in the absence 
of a magnetic background in the literature of the linear 
sigma model with quarks, coupled or not to Polyakov 
loops (e.g., Ref. [I3]). However, since this is a the- 
ory with spontaneous symmetry breaking (besides an ex- 
plicit breaking), a nonzero condensate modifies the quark 
mass, and vacuum subtractions are more subtle. At the 
end, there are finite logarithmic contributions that sur- 
vive renormalization, sometimes called zero-point correc- 
tions, but that are usually discarded phenomenologically 
in the LSMq. These contributions have been studied by 
the authors of Ref. [16] in the LSMq. Vacuum contribu- 
tions were considered at finite density in the perturbative 
massive Yukawa model with exact analytic results up to 
two loops in Ref. [33 IH] and, more specifically, in op- 
timized perturbation theory at finite temperature and 
chemical potential in Ref. |47j . also comparing to mean- 
field theory. More recently, this issue was discussed in a 
comparison with the Nambu- Jona-Lasinio model [691 . A 
very recent careful and detailed study in the (Polyakov 
loop extended) quark-meson model, with special atten- 
tion to the chiral limit, can be found in Refs. |48| . 

Thus, the one-loop correction from quarks to the free 
energy of the system is given by 



nqiB,T) 



r!p;-(s,T) . (31) 



f=u,d 



Below we discuss each term in this equation separately. 
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A. Vacuum contribution, Q 



vac 



The vacuum contribution can be expressed as the 
Heisenberg-Euler energy densit3Q 



log 



det(i^(«) - TO,) 



iVid I det{ip — rriq) 



(32) 



{qBf fds 



tanh s 



- 1 



-s7nl/{qB) 



where we have already accounted for the color degeneracy 
of the quark degrees of freedom. 

Since the covariant derivative p^'^^ contains not only 
the electromagnetic but also the non-Abelian gluon fields, 
one could expect that non-Abelian features should also 
appear in Eq. (32 1. Notice, however, that the Polyakov 
loop contribution is trivial since, for T = 0, it can be 
put away by a simple shift in the integral over the zeroth 
component of the momentum. So, for the vacuum con- 
tribution from quarks we can ignore its presence in the 
covariant derivative. On the other hand, as we discuss 
below, at non-zero temperature the presence of the gluon 
field leads effectively to a shift of the Matsubara frequen- 
cies, thus contributing to the paramagnetic (second) part 
of ([31]). 

The expression above corresponds, of course, to the 
quark vacuum correction in the presence of the external 
magnetic field 
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Nr. 



f=u,d 



(33) 



minus the vacuum correction in the absence of the field 



f=u,d 



where we have defined the integral 



(34) 



(35) 



^ ' I (27r)' 

and the quark mass in the presence of B 

Mlj^m)+2n\qf\B. 

One can control the divergences in the integrals above 
via dimensional regularization and compute the differ- 



(36) 



ence fi, 



{B) 



fig"-' in the MS scheme, obtaining 



Nr_ 
27r2 



E {qfBf[a-l,Xf) 

f=u,d 



^ In Ref. |68| one can find results for spinor and scalar determinants 
in the presence of external fields. 



1 

12 



{xj - a;/)log2;/ -f 
log 



2\qf\B 



(37) 



where Xf = mj: /{2\qf\B), is the derivative with re- 
spect to the first argument of the Riemann-Hurwitz C 
function |65| . and A is the renormalization subtraction 
scale. This expression is in agreement with the results 
of Refs. [49j [66], except for the scale-dependent term. 
Since the scale-dependent term is independent of the me- 
son fields and we are concerned only with the effective 
potential of the theory, we can discard it. Here our nor- 
malization is such that this correction vanishes in the 
limit of -B = 0. One arrives at the same result using the 
Heisenberg-Euler approach and imposing the normaliza- 
tion mentioned above. 

We can cast the correction to the effective potential 
in a dimensionless form normalizing every quantity by a 
mass scale. For this purpose, we choose v as our mass 
unit and define: 



b = 



eB 



t 



T 



(38) 



It is also convenient to write the quark electric charge 
as qf = rfcB sign(gy), so that = 2/3 and = 1/3. 
Then, the vacuum one-loop contribution to the effective 
potential is given by 



27r2 ^ 

f=u,d 

where the function F is defined as 



2rfb 



(39) 



F{x) = Ci-l,^f)-lix}-Xf) logxf + ^ , (40) 

and displayed in Figure [3] For small values of its argu- 
ment it behaves as 

F{x) = -0.165421 + ^x log x + 0{x) , (41) 

whereas for large a argument one finds 



Fix) = —(1 
^ ' 12^ 



log a;) + 0{x' 



(42) 



Since the constituent quark mass ruq is linear in cr, 
the external magnetic field modifies the classical poten- 
tial (17) for the a field. One can see from the behavior 



of this correction that its effect is to increase the value of 
the condensate and deepen the absolute minimum of the 
effective potential, as expected from previous results on 
magnetic catalysis. 

Finally, for the sake of completeness, we discuss the 
weak and strong field limits of the T = quark free 
energy ( 32 ) . The leading term in the weak- field limit 



{\qB\ <t:lr?y) is 



Nr 



3607r2 



qB 



1 + 



qB' 

m 2 



(43) 
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100 



FIG. 3. The function F{x), Eq. 1 40 1 



The strong- field hmit {\qB\ ^ mi) with on-shell renor- 
malization (A = m) is given by |68j : 



Qvac 



Nr. 



247r2 
N, 



{qBf In 



2 qBniq In 



167r2 



mt In 




(44) 



where Cq = e^/^-r = 2.51629, Ci e/27r = 0.432628, 
C2 = ^ 0.0505368, 7 = 0.577216 is the Euler con- 

stant, and A = 1.28243 is the Glaisher constani]^ 

Numerical tests show that only in the limit of very 
high fields, in units of v, the vacuum corrections bring a 
non-negligible modification to the effective potential. 



B. Paramagnetic contribution, 57^^"^^ 

The paramagnetic temperature-induced contribution 
comes from the second determinant ratio in Eq. ( 30 ) . In 



general, a quark determinant can be formally written as 
the product over the eigenvalues A„ of the corresponding 
Dirac operator, i.e.: 

det[ip + m,]T ^ n ^ ' (45) 

n 

where A„ satisfy the eigenvalue equation 



(46) 



It is important to stress that in Eq. ( 45 1 we omitted a T- 



independent constant which would anyway be subtracted 
at the end of the calculation. 



^ In Ref. 1521 . the quark vacuum term was not included, tlie 
cliarged pions vacuum contribution being considered instead. 
Both treatments increase the value of the condensate and deepen 
the effective potential. In this paper we choose to treat the 
mesonic fields classically. 



The algebraic equation for eigenvalues A„ reads as fol- 
lows: 



{Po 



-{2n + l~2s)\q\B ^0, (47) 



where q is the electric charge of the quark, n = 0, 1, 2, . . . 
is the quantum number that labels the Landau levels and 
s = ±1/2 is the projection of the spin of the quark eigen- 
state onto the z-axis. Here Aq is the ith diagonal com- 
ponent of the SU{3) gauge field. We assume no sum over 
color indices i = 1,2,3 unless explicitly indicated. From 



(471 one finds 



A^^=m, ± [ipo - Aff -pl-{2n + l- 2s)\q\B] (48) 

At zero magnetic field the contribution to the effective 
potential coming from (45) can be written as follows: 



lndet[i^ + m]T = 



(2^ 



In 



2^2 



, (49) 



where i labels the Matsubara frequencies uJi. The zeroth 
component of the momentum is related to the Matsubara 
frequency as follows: 



Po 



-ip4 



-iue, L0( = 2TrT{£+ 1/2) , (50) 



where the odd frequency takes into account the anti- 
periodicity of the quark fields along the temporal direc- 
tion. 

In the presence of a magnetic field B, the integral over 
the phase space in Eq. ( 49 1 is modified to 



d^p 
(2^ 



2tt ^ 

n=0 



dpz 
2tt 



(51) 



where the integer n labels the Landau levels 
Using the mapping 



49 



5l[) , the explicit expression for the 
and the eigenvalues ( 48 1 , we 



fermion determinant 
can rewrite the quark-induced potential (30) as follows: 

T 



^para 



(a,<f) 



In 



3d 

\q\BT 
2tt 



detT(«^^*^ - TO,) 
det(i^('3) - m,) 

3 00 

EE EE 



+00 



dpz 
27r 



In 



T 



^ln{Pz,(^) 
2^2 



(52) 



where the dispersion relation for quarks in the external 
magnetic field is 

w.n(Pz,a) = [m2(a)+p2 + (2n + l-25)|g|B]'/'(53) 
with constituent quark mass mq{a) — ga. 



With a suitable identification, 7 
LOsn/T, the sum over Matsubara frequencies in Eq. (52) 



A'l/T and 6 
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can be done explicitly with the help of the following for- 
mula (with arbitrary real 7 and 6): 

00 

^ ln[(27r£ + TT + + 6*2] =4^1n(27rj) (54) 

lei, j=i 

+9 + ln(l + e^'^e"^) + ln(l + e-'^^e-'^) . 

The first term in the right hand side of equation above 
is divergent. In the language of the effective poten- 
tial (52), the divergent term corresponds to a contribu- 



tion from massless fermions plus an additive constant. 
We ignore this term since it depends neither on chiral 
nor on confining variables. The second term leads to 



a T-independent contribution to ( 52 ) , which is already 



taken into account by the vacuum energy ( 32 1 . The last 



two terms lead to the regularized paramagnetic free en- 
ergy for the chiral and confining fields: 



+ C30 



2tt 



where 

Wt{<jJ,^) =Trc [ln(]l + $e- 



WT[c^s„(p.,a),*](,55) 



-I- c.c. 



(56) 



^Infl + 2e""/^ cos If ^ + e 



i=l 



-2w/T 



Here the trace Trc is taken in the color space and the 



angles tpi are defined in the diagonal representation ( 12 ) 



of the untraced Polyakov loop $ [the Polyakov loop is de- 
fined in Eq. ([I])]. Notice that the parameters (pi - which 
are the phases of the diagonal components of the un- 
traced Polyakov loop - enter the partition function ex- 
actly like imaginary chemical potentials. 



The paramagnetic contribution, given by Eqs. ( 55 ) and 



( 56 1 , can be simplified further by expanding the logarith- 

llowing 



mic function in Eq. (56) in the following series: 



In 1 + 2e 



COS if) 



k=l 

Next, we use the simple relation 

+00 

dpz exp ■ 



T 



T 



(57) 



(58) 



where Ki {x) is the modified Bessel function of the second 
kind (the MacDonald function) and order I [5S] . Finally, 
we express the free energy ( 55 ) in terms of the sums only: 



^para 



X^s„(cr) K 



± 1 n=0 fc=l 
k 



T 



Msn(o-) 



(59) 



where /i^n = tJsniPz = 0, tr) is the energy of the nth 
Landau level ( 53 1 at zero longitudinal momentum 



fisn{(r) - .9 + {2n + 1 - 2s)\q\B 



1/2 



(60) 



and the untraced Polyakov loop is defined in Eq. ( 12 ), so 
that 



Re [Tr*'^] =^cos(fc^O 



(61) 



The integer number k corresponds to the winding number 
of the Polyakov loops. 



TinL 



, 1.0 




IinL 



, 1.0 




FIG. 4. The sam e as in Fig. [2] but with the paramagnetic 
contribution ( 59 1 included. The potential is shown at two 



temperatures: T = 0.8 To (top) and T = 1.2 To (bottom). 
The magnetic field is given by \/ eB — 3T in both cases. 

It is also convenient to express the paramagnetic con- 
tribution to the effective potential in terms of dimension- 
less quantities as follows: 



27r2 

where we defined the dimensionless function 

00 3 

if(c, ,^2,6,0= E E EE 

/=ti,ds=±l/2n=0 i=l 



+2e-^^^+^^^^(«^cos(/.^) (63) 
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and the dimensionless version of Eq. (|60j), i.e. 

fisnf = [9^e + {'^n + l-2s)rfb]'^^ 



where we used Eq. (20) and defined to = Tq/v, besides 



(64) 



C. Paramagnetically-induced breaking of Z3 



One can see from Eq. ( 59 1 that the magnetic field dras- 



tically affects the potential for the Polyakov loop. Take, 
for example, the hmit of a very strong field, ^y\q\B 3> 
ruq = ga. Then, the leading contribution to Eq. ( 59 ) 
is given by the lowest Landau level, n = 0, with spins 
oriented along the field, s = +1/2, at lowest winding 
number of the Polyakov loop, k = 1. Thus, the leading 
term of the strong-field expansion of the paramagnetic 
potential (59) is 

,ga\q\BT^^ ( g_a\ 



Qpara 



ReL, 



(65) 



where the Polyakov loop variable is defined in Eq. ([T]). 
At fixed value of the magnetic field and temperature, the 
function (65) is a nonmonotonic function of the field a: 



it increases at small values of a and decreases at large 
values of a. This fact will be essential for our discussion 
of the role of the vacuum contribution introduced in the 
previous subsection: the vacuum contribution makes the 
expectation value of the a field larger, eventually dimin- 
ishing the role of the paramagnetic contribution ( 65 ) at 
large enough magnetic fields. 



Notice that the paramagnetic contribution ( 65 1 is not 



invariant under the Z3 symmetry, and therefore the para- 



magnetic contribution ( 65 1 deforms the potential for the 
Polyakov loops. So, it is clear that the magnetic field 
tends to break the Z3 symmetry and induce deconfine- 
ment in this model. 

In order to illustrate the effect of the paramagnetic 
term (59) we plot in Figure [4] the potential for the 
Polyakov loop (20) including the correction induced by 
the paramagnetic interaction ( |65[ ). In the figure. The 
magnetic field is \/eJ3 ^ 3T. We show the potential for 
two temperatures, T = 0.8 Tq ^J^d T = 1.2 To, in order 
to compare the plots with the pure gauge field potential 
displayed in Figure [2] The effect of the magnetic field 
is clearly seen in the deconfined phase: the paramagnetic 
interaction induces the expectation value of the Polyakov 
loop to be real-valued. 



D. Complete effective potential 

It is, of course, also convenient to define dimensionless 
expressions for the classical potential 



(66) 



where h = h/v^, and for and the potential for the 
Polyakov loops 

Vp{^i,<j>2,t,to) Vl{L,T) 

II = ' (67) 



expressing L and L* in terms of ipi and (j)2- 

The total (dimensionless) effective potential can then 
be expressed as 

Veffi^Al,^2,b,t) ^ V,i{0 Vp{<j,iA2.t,to) 

y4 y4 y4 

. Kac(e,fe) , l^P^'-"(C,^l,c/.2,6,0 

(68) 

and it is this form of the effective potential that we use 
to determine the phase structure of the PLSM^ effective 
field theory in the following section. 



IV. PHASE STRUCTURE 

In our numerical analysis we have considered the fol- 
lowing situations: 

(A) The linear sigma model with thermal corrections, 
without the coupling to the Polyakov loop. This is 
a standard analysis and has been performed previ- 
ously in 9 - 20 , and it is shown here for the sake of 
completeness; 

(B) The linear sigma model at finite temperature cou- 
pled to the Polyakov loop. At this point the effects 
of the magnetic field were not taken into account 
yet and we investigated only the effects of the in- 
teraction between the two order parameters; 

(C) The linear sigma model in a magnetic background 
at zero temperature. In this case the Polyakov loop 
is excluded by construction. The magnetic field 
manifests itself through vacuum corrections from 
the quarks; 

(D) The complete scenario: the linear sigma model cou- 
pled to the Polyakov loop with thermal corrections 
in the presence of a magnetic background. 

In the latter case we also investigated the effects of in- 
cluding vacuum corrections or not. The discussion about 
the inclusion or not of such terms is still inconclusive, 
not being clear which treatment is closer to the original 
theory, QCD. We obtained phase diagrams in the (T, B) 
plane for both cases and verified that the inclusion or 
not of vacuum corrections is crucial to the structure of 
the phases, influencing not only the nature of the transi- 
tions, as it was first observed in 69J, but also changing 
the qualitative behavior of the transition lines. 



A. B = 0,T / 0, (pi,2 = 

This case is described by the linear sigma model cou 
pled to quarks. We take into account thermal correc 
tions and disregard coupling to the Polyakov loop (12 1 
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The chiral condensate is the approximate order parame- 
ter for the chiral transition, and for low temperatures it 
is nonzero, indicating symmetry breaking. The value in 
the vacuum is adjusted to correspond to the pion decay 
constant, /^r, as described in section II C As the tem- 



perature is raised the minimum approaches zero, and for 
high enough T the symmetry is recovered. As the po- 
sition of the minimum as a fonction of temperature is 
continuous, and so is its derivative, the transition is a 
crossover. 




FIG. 5. Effective potential of the linear sigma model with 
thermal corrections for several values of temperature. The 
transition is a crossover. 



B. B = 0, T / 0, V3i,2 / 

In this subsection we analyze the interaction between 
the chiral condensate and the Polyakov loop, without the 
presence of the magnetic field. The treat the model at fi- 
nite temperature because of the presence of the Polyakov 
loop ( 12 ) in the formulation of the model. The expression 



for the thermal effective potential is similar to Eq.(55l 



q 



dpi 
(27r)3 



W^s=o[w(Pz,fT),$] 



(69) 



but 



where Wb=o has an expression similar to Eq.(56), 
the frequency u has no dependence on the magnetic field, 

and is given hy oj — \Jp'^ + 

Due to the interaction, the Z3 degeneracy is broken 
and the minimum along the real axis becomes the global 
one. Both transitions, chiral and deconfinement, are con- 
sistent with a crossover-type transition since the order 
parameters behave as smooth functions of the temper- 
ature. Fig. [6j The slopes of both order parameters be- 
comes steepest at the same temperature Tc ~ 215 MeV, 
clearly showing the tight relation between the order pa- 
rameters. 




200 250 300 

T (MeV) 



400 



FIG. 6. Expectation values of the order parameters for 
the chiral and deconfinement transitions as functions of the 
temperature. The filled circles represent the a-condensate, 
and the empty circles stand for the expectation value of the 
Polyakov loop. Both lines are smooth functions of the tem- 
perature indicating the presence of a crossover-type smooth 
transition. 



C. r = o 

At zero temperature the Polyakov loop is absent by 
construction. The only pieces of effective potential that 
contribute are the linear sigma model potential and the 
vacuum correction from quarks. The magnetic field in- 
fiuencing the system through the latter. We plot in Fig. [7] 
the effective potential for different values of the external 
magnetic field. We can see that as the field increases in 
magnitude the potential becomes deeper and the value 
of the condensate raises, in accordance with results at 
zero temperature that indicate an enhancement of chi- 
ral symmetry breaking. This effect - known as magnetic 
catalysis - is in accordance with previous result from two 
of us [52]. In Fig. [7] we subtracted the value of the poten- 
tial with ^ = 0, in order to have a better comparison for 
different values of magnetic field, and AV is the resultant 
potential. This procedure will be applied in all plots of 
the effective potential along this section. 



D. 



T/O, 0/0 



Finally we analyze the full potential. It contains the 
linear sigma model, the pure gauge potential, thermal 
corrections including interaction between the two fields, 
and vacuum corrections from the quarks. The effects of 
including or not the vacuum terms were also considered 
and have shown to be crucial for the structure of phases. 
The relevant potential is given in Eq. (68). This potential 
is a function of the parameters cr, 931, and (^2, where we 
used the parameterization (12 1 and set Lp^ = —^pi — ip2- 
As shown in Fig. [2] for pure gauge the potential - written 
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Without vacuum corrections 



100 



FIG. 7. Zero-temperature efTective potential (68 1 with vac- 
uum corrections included is shown as a function of the chiral 
variable = a/v for different values of external magnetic field. 



in terms of L and L* - has three degenerated minima. It 
is known that in the presence of quarks this degeneracy is 
broken and the global minimum will be in the direction of 
Re[L]. Since we are interested only in the minima of the 
potential, we can simplify our analysis by setting Im[L] 
to be zero. 



Without vacuum corrections 



10 



S 
o 

It 
m 



> 



5 - 











— T=150MeV 






- T=185MeV 






■ ■ T=205MeV 






, . • ' 


1 

1 

/ 




/ y 
/ / 

^ / 








1 




1 



0.5 



1.5 



FIG. 8. Complete effective potential without vacuum correc- 
tions from quarks as a function of ^ = ct/h. The value of <f 
was chosen to take the minimum for each temperature. The 
barrier between the two minima indicates a first-order chiral 
transition. 

Let us consider first the potential without vacuum cor- 
rections. We evaluate it for different values of magnetic 
field and temperature. For each set of values of T and 
B we minimize the potential with respect to a and 
Then we use this value of to plot the effective po- 
tential as a function of a only. Figure [8] illustrated the 
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FIG. 9. Complete effective potential without vacuum correc- 
tions from quarks as a function of Re\L\. The value of a was 
chosen to be the minimum for each temperature. The differ- 
ent curves are then slices with different cr, representing a first 
order transition, even if there is no apparent barrier in the 
plot (the first order jump in the chiral condensate makes the 
expectation value of the Polyakov loop discontinuous as well) . 



evolution of the chiral transition at the magnetic field 
eB — 10m J. Initially the condensate has a non-zero 
expectation value and this value decreases as the tem- 
perature increases. Then, a local minimum appears for 
cr = 0, and above a certain critical temperature it be- 
comes the global minimum. The barrier between the two 
minima causes a discontinuity in the expectation value 
of (T, pointing out to a first order character of the chi- 
ral transition, in accordance with [52). In the same way 
as in the case i? 7^ 0, T = 0, the potential becomes 
deeper for high values of the magnetic field. However, 
the thermal corrections, responsible for restoring chiral 
symmetry, also acquire a higher magnitude, bringing the 
chiral critical temperature down. 

InFig.|9] we fix CT at its minimum for a set of values of T 
and i?, and plot the potential as a function of We can 
then follow the evaluation of the expectation value of the 
parameter Ke\V\. At low temperatures it is zero, indicat- 
ing confinement, and above a critical temperature it goes 
asymptotically to one. At the critical value of the tem- 
perature, where the chiral transition occurs, the variable 
(L) also presents a discontinuity, so that the first-order 
chiral transition is accompanied by a first-order decon- 
finement transition for large strengths of the magnetic 
field. Thus, if the vacuum corrections are not taken into 
account, then the magnetic field enhances the order of 
the phase transition and the transition changes from the 
crossover type (realized at £? = 0, Section IV B I to the 
first order type. 



The phase diagram for the potential described above is 
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Without vacuum corrections 



200 - 



■ Chiral transition 
Deconfinement transition 




125 - 



FIG. 10. Phase diagram for the complete potential without 
the vacuum corrections coming from quarks. The chiral (black 
line) and deconfinement (light red line) transitions happen at 
the same temperature, leading to a common line. Both critical 
temperatures decrease with the increase of B. The lines mark 
first-order transitions, except for very small values of B, when 
the transitions become crossovers. The exact value of field for 
which it happens is not clear from our numerical analysis. 



shown in Fig. 10 The critical temperature for the chiral 
and deconfinement transitions are the same for different 
values of magnetic field, resulting in two phases: a con- 
fined phase with broken chiral symmetry, and a decon- 
fined phase with restored chiral symmetry. The common 
critical line of these phase transitions is of first order for 
all values of the magnetic field B except for very small 
values of B. As the magnetic field decreases, the transi- 
tion turns back into a crossover reproducing the physical 
situation with B = that was discussed around Fig. [6j 

The inclusion of vacuum corrections from quarks 
changes the phase diagram dramatically. Following the 
same procedure describe above, we plot the potential as 
a function of a in Fig. 11 and as a function of Re[L\ 



in Fig. |12| The corresponding phase diagram is shown 



in Fig. [13) The vertical green line indicates the value of 
magnetic field that is expected to be reached in heavy-ion 
collisions in the ALICE experiment at the LHC [7 . 

The presence of the quark vacuum term changes the 
order of the chiral transition, turning it into a crossover. 
In the case without the vacuum correction, the transi- 



^ The starting point of the curve — given by the critical temper- 
ature at zero chemical potential - is somewhat higher compared 
to the value obtained in recent numerical calculations of lattice 
QCD with two flavors of light fermions (see, for instance, the 
review 1701 V Generally, effective models like our PLSMg or the 
PNJL model give approximate answers for particular quantities 
while being able to predict a general picture at a very good qual- 
itative level. 
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FIG. 11. Complete effective potential with vacuum correc- 



tions (68 1 coming from quarks as a function of ^ = cr/u. The 



transition is a crossover. 



With vacuum corrections 
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FIG. 12. Complete effective potential with vacuum correc- 
tions from quarks as a function of i?e[L] for B = 20m^ and 
several temperatures. 



tion was of first order due to thermal effects. The ther- 
mal contribution generates a minimum at ^ = being 
rather sensitive to the value of the magnetic field. In the 
present case, however, the vacuum term has a magni- 
tude that is larger compared to the thermal piece. Thus, 
the vacuum contribution dominates the position of the 
minimum, and it is the vacuum part that is responsi- 
ble for driving the transition to a crossover. Moreover, 
the chiral-violating property of the vacuum quark contri- 
bution overwhelms the chiral-restoration tendency of the 
thermal part so that the critical temperature of the chiral 
transition - calculated by taking into account both vac- 
uum and thermal effects - becomes an increasing function 
of the external magnetic field. 

It is interesting to recall that the thermal contribu- 
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FIG. 13. Phase diagram for the complete potential with the 
vacuum corrections from quarks. The critical temperatures 
of the deconfinement (the dash-dotted line) and chiral (the 
dashed line) transition coincide at B = and split at higher 
values of the magnetic field. A deconfined phase with broken 
chiral symmetry appears. The green solid vertical line is the 
magnitude of the magnetic field that expected to be realized 
at heavy-ion collisions at the LHC [7]. 



increases; and that (ii) at some critical value of the mag- 
netic field, eB ~ 24171^, color confinement may disappear 
as illustrated in Fig. [T] On the contrary, we show that 
the inclusion of the mentioned vacuum corrections make 
the critical deconfining temperature to be an increasing 
function of the magnitude of the magnetic field up to 
eB ~ 50to^. We have not found a signature of a criti- 
cal magnetic field that would lead to the magnetic-field- 
induced deconfining transition in either scenario, with 
and without vacuum corrections. 



V. CONCLUSION 

In this paper we studied the influence of a strong mag- 
netic field background on confining and chiral properties 
of QCD, using the linear sigma model coupled simultane- 
ously to quarks and to the Polyakov loop (PLSMg). This 
model associates electromagnetic, chiral and the confin- 
ing degrees of freedom in a natural way. 

In the confining sector, the strong magnetic field af- 
fects the expectation value of the Polyakov loop, which 
is an approximate order parameter for the confinement- 
deconfinement phase transition at finite temperature. 
We found that the contribution from the quarks to the 
Polyakov-loop potential has three features: 



tion of quarks tends to destroy the confinement at strong 
magnetic fields. Fig. 
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contrary to the confinement- 
enhancing property of the vacuum contribution. This 
effect is due to the fact that the vacuum contribution 
makes the value of the field a larger (as one can see in 
Fig.[7|), suppressing the Z3 violating paramagnetic contri- 



bution (65 1 at large enough magnetic fields. Thus, large 



magnetic fields tend to enhance the confinement property 
of the vacuum. As in the case of the chiral transition, the 
vacuum contribution exceeds its thermal counterpart so 
that the deconfinement transition temperature appears 
to be a weakly increasing function of the magnetic field. 
Fig. [TSj The confinement transition is also a crossover. 

The lines of the confinement and chiral transitions co- 
incide for small values of B, and go apart as the field in- 
creases in magnitude. Results from lattice QCD indicate 
that in the absence of the magnetic field the deconfine- 
ment transition and chiral symmetry restoration happen 
in the same narrow temperature interval |71j . Our cal- 
culations indicate that the presence of a strong magnetic 
field should inevitably split these transitions in this case. 

The deconfinement transition occurs at lower temper- 
atures compared to the temperature of chiral restora- 
tion. Thus, in a strong magnetic field a deconfined 
phase with broken chiral symmetry should appear. The 
splitting between these transitions becomes wider with 
the increase of the magnetic field B. Some other ap- 
proaches IH] - that were not taking into account vac- 
uum corrections - demonstrate that (i) the deconfine- 
ment temperature drops down as the magnetic field B 



1. the presence of the magnetic field breaks the global 
Z3 symmetry and makes the Polyakov loop real- 
valued (this effect is seen in the Polyakov loop po- 
tential, Fig.|4]); 

2. the thermal contribution from quarks tends to de- 
stroy the confinement phase by increasing the ex- 
pectation value of the Polyakov loop; 

3. on the contrary, the vacuum quark contribution 
tends to restore the confining phase by lowering 
the expectation value of the Polyakov loop. 

The vacuum correction from quarks has a crucial im- 
pact on the phase structure. If one ignores the vacuum 
contribution from the quarks, then one finds that the 
confinement and chiral phase transition lines coincide. 
Fig. 
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and in this case the increasing magnetic field 
lowers the common chiral-confinement transition temper- 
ature. However, if one includes the vacuum contribution, 
then the confinement and chiral transition lines split, and 
both chiral and deconfining critical temperatures become 
increasing functions of the magnetic field. Fig. [l3| The 
vacuum contribution from the quarks affect drastically 
the chiral sector as well. 

Our calculations also show that the vacuum contri- 
bution seems to soften the order of the phase transition: 
the first-order phase transition - which would be realized 
in the absence of the vacuum contribution - becomes a 
smooth crossover in the system with vacuum quark loops 
included. 
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It is important to stress that in a strong magnetic field 
a deconfined phase with broken chiral symmetry will ap- 
pear in the scenario with quark vacuum corrections. The 
splitting of the confinement and chiral transitions - as 
shown in Fig. 13 - may be substantial at a steady mag- 



netic field of the magnitude that is expected to be realized 
in heavy-ion collisions. For example, in the ALICE ex- 
periment at the LHC facility, the magnetic field may peak 
around the value eB ~ 15m^ 0. We find (see Fig. 13) 



that the splitting between the critical temperatures of 
the confinement and chiral crossovers in a constant field 
of the typical LHC magnitude may reach the noticeable 
value 



contributions |72) . 

Note added. After this paper was submitted to the 
journal, a new preprint on a similar subject has ap- 
peared [73]. The authors have found a splitting of the 
deconfinement and chiral transitions in the magnetic field 
background in the framework of a Polyakov-loop ex- 
tended Nambu-Jona Lasinio model with eight-quark in- 
teractions. This independent result is an agreement with 
our conclusion on splitting of these transitions in the lin- 
ear cr-model (with the vacuum corrections included). 
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